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Abstract: We present a detailed study of the Entanglement Entropy for the confin¬ 
ing Klebanov-Strassler background coupled to a large number of dynamical flavors in the 
Veneziano limit. As we vary the number of the massless flavors the behavior of the entropy 
strongly depends on the way we fix the integration constant of the warp factor, that is 
related to the glueball scale. In the case of massive flavors, the mass of the flavor branes 
introduces another scale in the background and the entropy undergoes two first order phase 
transitions. The competition between the glueball and the quark scales will lead to a crit¬ 
ical point where one of the phase transitions degenerates to a second order one. We have 
calculated the critical exponents and have found that they are independent of the number 
of flavors and different from the mean filed theory expectations. 
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1 Introduction 

One of the most inspiring conjectures in modern theoretical physics is the AdS/CFT corre¬ 
spondence [1] (for a set of pedagogical introductions see [2-4]) which, in its original form, 
claims the equivalence between the maximally supersymetric field theory in 4 dimensions 
and type-IIB superstring theory on AdS§ x S 5 . The original AdS/CFT correspondence can 
be generalized to include field theories with less supersymmetry, as well as field theories 
living in a different number of spacetime dimensions. The main interest in this kind of 
dualities is due to the fact that they provide a means to probe the dynamics of strongly 
coupled gauge theories. This is so because the duality is between strongly coupled gauge 
theory, on one hand, and weakly coupled string theory, on the other and vice-versa 1 . 

More precisely, one of the long-standing problems has been that of describing the 
confinement in theories like QCD. The aim of this work is to clarify how the presence 
of dynamical quarks in a theory affect the phenomenon of confinement. To this end, we 
choose to study the wrapped deformed conifold background (KS model [7]) both with 
massless and massive backreacting flavors ([8] & [9]). The KS background corresponds to 
the low energy theory of N regular and M fractional D3-branes on the deformed conifold 
Z\Z 2 — Z 3 Z 4 = e 2 . The dual field theory is a N = 1 4-dimensional gauge theory with initial 
gauge group SU(N+M) x SU(N), where N = nM and n is an integer. The bifundamental 
matter fields A, B interact through the superpotential W = e kl AiB^AjBi and transform 
as doublets under the SU (2) x SU (2) isometry group of the deformed conifold. This theory 
develops a Seiberg duality cascade which terminates after n — 1 steps, as soon as the gauge 
group has been reduced to SU(2M) x SU(M). The important point is that the KS theory 
exhibits confinement due to the formation of a gluino condensate < AA > A QCD ' I* 

lr This fact is at the same time a disadvantage since it is very hard to prove the conjecture. In the case 
of the original AdS/CFT scenario progress on the comparison of the spectra of the two theories has been 
possible due to integrability [5]. More recently, there was also some progress on the comparison of higher 
point correlators, see [6] and references therein 
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apparent that the confining scale A qcd is related to the deformation parameter e which 
breaks the scale invariance of the supegravity solution. 

The addition of fundamental degrees of freedom to the theory in the Veneziano limit 
can be achieved through the smearing of Nf flavor D7-branes along the internal directions 
(see [10] for a review on the smearing technique and [11] for some more smeared solutions). 
This construction introduces a second scale to the theory related to the bare mass of 
the quarks. It is the competition between these two different scales which will give the 
interesting phenomenon of two quantum first order phase transitions in the connected 
part of the Entanglement Entropy (EE). One of the phase transition is associated to the 
confinement-deconfinement phase transition and is present in the unflavored KS theory, 
while the other is related to the presence of dynamical quarks. As we decrease the mass 
of the quarks, approaching the glueball scale, we reach to a point where one of the phase 
transitions degenerates to a second order one. Beyond this point we end up with a single 
first order phase transition. For the details of the construction of the supergravity solutions 
we refer the reader to section 3, as well as to the original literature [7-9]. 

From what is said above it is inferred that the observable that we will use as an 
order parameter for the confinement-deconfinement phase transition will be that of EE. 
In particular, we will calculate the EE between region A and region B in the strongly 
coupled KS background with Nf dynamical flavors. More precisely, we will divide one of 
the spatial directions into a line segment of length L, which will call region A, and its 
complement, which will call region B. Then the EE between region A and region B can be 
used as a measure of the effective degrees of freedom at the energy scale 1/L. In [12-14] 
a holographic prescription for the computation of the EE was proposed, which consists of 
finding a minimal surface with a boundary that coincides with that of region A. Actually 
there are two different minimal surfaces with the same end points at the boundary. One 
is with two disconnected straight lines extending from each of the boundaries to inside the 
bulk while the other is a connected surface that connects the two end points of region A. 
In a generic confining theory, there is a critical length L c above which the disconnected 
lines are favored, while below it the connected curve that is favored. Thus in a large N 
confining theory the behavior of the EE becomes trivial in the IR (L > L c ) in the sense 
that 8 lS ~ N° << N 2 while it is of order N 2 for L < L c . This change of behaviour can 
be thought as the confinement-deconfinement phase transition. 

The plan of the paper is as follows. In Section 2 we assemble all the necessary tools 
for the holographic computation of the EE. In Section 3 we discuss the construction of 
the supergravity solutions corresponding to the inclusion of massless & massive dynamical 
flavors to the KS background. In Section 4, we analyze the behavior of the length of the 
strip and of the EE as we vary the number of the massless flavors. Depending on which 
scale we keep fixed we will find different qualitative dependence of the afore-mentioned 
observables on the number of flavors. 

In Section 5, we perform a similar analysis for the massive flavors. The EE undergoes 
two first order phase transitions, one associated with the confinement-deconfinement phase 
transition while the other with the presence of the dynamical flavors. Below some critical 
value of the quark mass one of the first order phase transitions degenerates to a second 
order one. Calculating the critical exponents we find that their values are independent of 
the number of flavors and different from the mean-field theory expectations. Finally, in 
Section 6 we draw our conclusions. 
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2 Entanglement entropy computation in a confining background 


In this section we will provide a short review of all the necessary results from [14] (see 
also [15]), in which the authors generalized the Ryu-Takayanagi conjecture [12, 13] for the 
gravity duals of confining large N c gauge theories. 

In a quantum field theory the EE between two complementary spatial regions A & B 
is the entropy measured by an observer in region A, who has no access to the degrees of 
freedom in region B. In [12-14] a holographic computation was proposed, which consists of 
finding a minimal surface with a boundary that coincides with that of region A. Then, the 
EE between region A and the complementary region B is given by the following expression 

s = i (2 - 1) 

where Gqo is the ten-dimensional Newton constant & Gs is the induced string frame metric 
on E. The calculation of the EE is realized through the minimization of the action (2.1) 
over all the surfaces that approach the boundary of the entangling surface. 

A generic expression of the string frame metric for the type of backgrounds that we 
will consider is 

ds 2 = a(r) [^(r)dr 2 + dx^dx^] + gij d6 l dO^ , (2.2) 

where x^ (g = 0,1,..., d) parameterize M rf+1 , (i = d + 2,..., 9) are the 8 — d internal 
directions and r is the holographic radial coordinate 2 


ta < r < r* . 


(2.3) 


The background is supplemented by a dilaton (that we will denote by 4?) and some fluxes 
(which will not enter in the EE calculation). The volume of the internal manifold is 


3 -d, 


Vint — 


Q dd ' l y de t[d 


i —1 


(2.4) 


In [14] the authors focused their attention on entangling surfaces that are strips of 
length L. Analyzing the equation of motion they found that it has two independent so¬ 
lutions corresponding to two local minima of the action (2.1). The first is a disconnected 
surface, consisting of two cigars that extend in and separate in the remaining di¬ 

rection by a distance L. The second is a connected surface, where the two cigars are 
connected with a tube whose width depends on L. The length of the line segment for the 
connected solution is 


l(t 0 ) = 2 r 

J To 



(2.5) 


where tq is the minimal value of r along the connected surface in the bulk and we have 
also defined the following quantity 3 


H(t) = e~^V? nt a d . (2.6) 

2 ta is zero for the KS geometry while there is a finite UV cut-off (due to the presence of a Landau pole) 
when we consider the KS background with massless/massive dynamical flavors. 

li The complicated integral in (2.5) can be approximated unexpectedly well with a combination of the 
background functions that appears in equation (3.1) of [15]. This claim has been checked for some of the 
results that we will present in the rest of this paper. 
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The EE of the connected surface as a function of tq is given by 


Sc{tq) = [ dT 

■0 J To 


2 Gio 


\ 


1 _ H(tq) 
1 H(r) 


(2.7) 


while the EE of the disconnected surface is independent of tq and given by 


Sd(t 0 ) = 


Vd~i r 

2 Gio Jt a 


dr \//3(r) H(t ). 


( 2 . 8 ) 


Even if the EE is in general UV divergent, the difference between (2.7) and (2.8) is finite 


2 Guo 
Ki-i 


5(r 0 ) 


2 Gio 


( 5 C - s D ) = r 

J TQ 


dr 


N 


1 _ ^(a) 


dr y/p{T)H{r), 


' TA 


(2.9) 


and it is this quantity that we will calculate in the specific case of the KS background with 
massless/massive dynamical flavors. 


3 Klebanov—Strassler background with dynamical flavors 

In this section we will review the two known solutions of [8] and [9, 16] for the addition 
of massless and massive flavors to the Klebanov-Strasser background [7]. We will arrange 
the constants in such a way that the limit from the backreacted to the unflavor solution is 
smooth, something that will be useful when we will perform the calculation of the EE. 

3.1 Addition of massless dynamical flavors 

The background we will construct accounts for the addition of fully backreacted (fractional 
and regular) D3-branes and (smeared flavor) D7-branes on the deformed conifold. The 
action consists of the usual color part (type IIB supergravity) and the flavor part (Dirac- 
Born-Infeld & Wess-Zumino terms). The ansatz for the metric (in the Einstein frame) we 
will adopt is the following [8, 9] 

ds 2 = /i _1//2 (r) dx^ dx^ + /r 1//2 (r)dsg with (3.1) 

= I e2G> ^ +* 2 > +e2Gl (' - db) W' + «+(‘ + db) « + <&) • 

where 

91 = 72 ^ 2 “ CT2 ) ’ 93 = 72 ^ 2 + a2 ) > 

92 = 72 (~ Ul + (Tl ) ’ 54 = ^(^i + ai) & g 5 = u; 3 + a 3 , (3.2) 

with 

ui\ = sin if; sin 62 dp 2 + cos 1/2 dO 2 , 0 J 2 = — cos ^ sin 62 dp >2 + sin 7 d 62 , 

oo 3 = d'l/j + cos 62 d ! ^2 , cr 1 = d 9 \ , 02 = sin^id^i & a 3 = cos 61 dipi. ( 3 . 3 ) 
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The range of the angles is G [0, 47t), ipi G [0, 27t), 0, G [0,7r] and t G [0,oo). The 
expressions for the dilaton and the (RR and NS) forms are (in units of g s = 1 Sz a! = 1) 


Fr, = dh 1 (r) A dx u A • • • A dx 3 + Hodge dual, = 4>(r), 


B 2 = 


M 

2 V 

H 3 = — 

2 

M 


fg 1 Ac / 2 + kg 3 A g 4 


F 1 = ^ g 5 
47r 

1 


dr A (/ g 1 A g 2 + k g 3 A g 4 ) + -{k - f) g 5 A ( 5 1 A £ 3 + g 2 A g 4 ) 


= y{^ A V F+ ^ f ^ gl ^ 9 2 + {l~ F +^k)g 3 Ag 4 \ + Fdr A (g 1 A 5 3 + g 2 A 5 4 )} , 


where M is the fractional D3-brane Page charge. The Gi = G,(r) (i = 1,2,3), h = h(r), 
f = /(t), k = k(r) and F = F(t) are unknown functions of the radial coordinate, while 
the dot denotes derivative with respect to r. 

The Bianchi identities for H 3 and F 3 are automatically satisfied by the proposed ansatz 
while the one for T 5 can be reduced to the following first order differential equation 


4 he 2Gl+2G2 


— M 


f-(f-k)F + Ffk 


(3.4) 


Applying the fermionic supersymmetric variations we obtain a system of BPS equations 
that fully determine the unknown functions of the background. The BPS equations for the 
functions of the metric are 


G\ 

G2 

G 3 

and for the dilaton it is 


1 2G 3 —Gi- 

9 


(3.5) 


N 


$ = — L 


f '!> 

4 e 

47T 




Nf , 

= 1 + 4 ^ (T * " T) 5 


(3.6) 


where we have introduced an extra scale r*, at which we fix the value of the dilaton. Notice 

that the solution in (3.6) is defined up to r < t^p, where tip = r* + is the point 

that the dilaton blows up. In order to write in a compact form the solution of the BPS 
system in (3.5), we introduce the following set of auxiliary functions A(r) & 7 /(r ) 4 


A (r) = 


1 


2 1 / 3 sinhr - 


sinh 2 r — 2 r + r/(r) 


■ , d 

with — 
dr 


-<s> 


N 


rj) = -p- (sinh 2 r — 2 r) 


(3.8) 


4 We did not integrate the differential equation for g in (3.8) in order to make smooth contact with the 

massive case that we will present in the next subsection. If we integrate and substitute the expression for 
77 in the one for A we have 


A(t) 


2 1 / 3 sinh' 


sinh 2 r — 2 r + 


Nf cosh 2 r — 2 r 2 
1 + Tp (r*~ t) 


87 r 


(3.7) 
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The metric functions Gi in terms of the expression for the dilaton (3.6) are 


e 


2Gi 


4 


£3 $ 



sinh 2 r . . . 
coslG 4 T ’ 


„2G 2 


€ 3 i * / \ o3 4 6 3 

-e 3 coshrA(r), e 3 = - e5 , 


(3.9) 

where e is an integration constant. Supersymmetry imposes the following constraints on 
the functions that appear in the expressions for the fluxes 


k = e* (f + ^ /) coth 2 T - , / = e' 1 ’ (l — F + ^ fc) tanh 2 ^ , F = ^ e $ (/c - /), 

(3.10) 

which are solved as follows 


f = e 


t coth r — 1 
2 sinh r 


(cosh r — 1 ), k = 


t coth r — 1 
2 sinh r 


(coshr + 1), F = 


sinh t — t 
2 sinh r 


(3.11) 

where the dependence on the number of flavors is through the equation for the dilaton 
(3.6). Combining (3.9) with (3.11) and (3.4) we can immediately solve for the warp factor 
h(r). The integration constant can be fixed in different ways and in section 4 we will see 
the effect of two alternative possibilities on the qualitative behavior of the EE. 


3.2 Addition of massive dynamical flavors 

The BPS equations between the massless [ 8 ] and the massive [9] case are unchanged, up 
to the substitution 

Nf => Nf(r). (3-12) 

This is the case since the supersymmetric variations contain only expressions of the F t ‘s 
and not of their derivatives, despite the difference between the modified Bianchi identities 
of the massive and the massless case. The extra ingredient in order to solve the BPS 
equations is to determine the distribution function of the D7-brane charge density by a 
microscopic calculation, which was performed in [9, 16]. The outcome of that analysis is a 
very complicated expression that looks like a smoothed-out Heaviside step function 

N f (r) = N f 0(r - r q ), (3.13) 

where T q is related to the dynamical quark mass 5 . Comparison between calculations using 
the exact expression of the brane distribution function and the Heaviside approximation 
(3.13), of the quark antiquark potential [20] and of the meson spectrum [21] in the case 
of the massive Klebanov-Strassler and Klebanov-Witten backgrounds, have shown that 
qualitatively the physics remains unchanged. Based on this observation we will perform 
the calculations of the EE using the heaviside approximation, since additionally using the 
exact expression of the brane distribution function makes the computation extremely time 
expensive. 

The smeared set of massive flavor branes produces a spherical cavity that effectively 
splits the solution into two regions. In the inner region, where r e [0,r g ], the effective 
charge for the D7-branes is zero and imposing regularity at the level of the cavity the 
solution is a mild modification of the known KS unflavored result. In the outer region, 

5 At the probe level the radius of the spherical cavity, that the flavor branes distribute, correspond to 
the bare mass. When backreaction is taken into account see [17-19]. 
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where r > r q , we have to solve the new system of BPS equations with the substitution 
(3.12). 

Starting with the dilaton equation (3.6) we can easily integrate and, after imposing 
regularity at the cavity, obtain the following expression 

1 + £ (r* - T g ) if r <r q 
1 + £ ( T * - r) if r > r q 

Following the same line as in the massless case we are able to solve the massive BPS system 
and the only function that needs to be determined is t](t). From (3.8) we can see that 
this function has a constant value inside the cavity, where the unflavored solution lives. 
Imposing continuity at r = r q we can fix this constant to zero, so the expression for r/(r) 
becomes (in the Heaviside approximation) 



(3.14) 


r/(r) 


0 

Nf cosh 2t—2t 2 — cosh 2t 9 H-2t 2 
87T 1 , N f / \ 

1+ 4lr ( T * -T ) 


if r < Tq 
if T > T q . 


(3.15) 


Taking into account that the expressions for the fluxes are still given by (3.11), having 
determined the dilaton and the function r), from (3.14) & (3.15) respectively, completely 
solves the problem. As in the massless case the expression for the warp factor h(r) is 
obtained by integrating (3.4) and imposing continuity at the level of the cavity. 

In order to check the regime of validity of the supergravity dual of the current section, 
we can study, in the string frame metric, the Ricci scalar R and the square of the Ricci 
tensor RmnR AIN ■ It is not difficult to check that 


$ 1 o , M 2 

e ~ —— <fc n ~ —— 

N f N f 


R 


N f 


yfh^ M 


(3.16) 


and therefore the solution is reliable in the interval 1 -C Nf <C M. This calculation is valid 
for both the massless & the massive solution. 


4 Entanglement entropy massless flavors 

In this section we will place a strip of length L at a finite UV cut off scale of the Klebanov- 
Strassler geometry with massless dynamical flavors and calculate the EE difference, as it 
is defined in (2.9). Since we want to emphasize the effect of the dynamical flavors we fix 
a comparison scheme in which we vary the number of flavors, while we will keep fixed 
the UV scale r* and the integration constant e (defined in (3.6) & (3.9) respectively). An 
alternative choice would be to keep fixed the Landau pole, as it is defined below (3.6), 
but to ensure a smooth zero flavor limit (unflavored solution) we need a flavor dependent 
Landau pole 6 . 

As we mentioned in subsection 3.1 another constant that needs to be fixed is the 
integration constant that comes from the equation of the warp factor h(r), namely ho. 
This constant determines the scale for the glueball and KK masses, since m g i ue 1/ y/ho- 
We will consider two specific choices of fixing the constant ho ‘: 

6 For the possible physical implications of the different comparisons schemes see [22] . 

' In [9] the authors examined the effect of two alternative possibilities on the Wilson loops calculation. 
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• Case I. We fix ho, in analogy with the unflavored KS solution [7], by requiring that 
h(r) goes to zero at some finite UV cut off (before reaching the Landau pole that the 
supergravity approximation is no longer valid) for any number of flavors 8 . In this case 
ho(Nf), as can be seen from the left plot of figure 1, is a monotonically decreasing 
function and as a consequence the glueball mass increases with the number of flavors. 

• Case II. We fix ho by imposing that the value of h(r) at the origin is the same for 
any value of Nf, see the right plot of figure 1. Here we invert the strategy that we 
followed in Case I, and since the numerical value of the glueball mass may depend 
on other parameters of the geometry besides ho, we will obtain a glueball mass that 
decreases with the number of flavors. The price we have to pay for such a choice is 
that the warp factor is finite in the UV, which amounts to deforming the QFT dual 
with the insertion of an irrelevant operator (see the discussion in [23]). 


h.icr 4 h.icr 4 



Figure 1. We have plotted the warp factor as a function of the holographic coordinate for the two 
different choices of fixing the integration constant ho. The black curve is for zero flavors and in 
both cases matches the computation of [14]. The red(blue) curve is for Nf = 0.1 (Nf = 0.5). In the 
numerical evaluation we have used e = 1 , M = 10 8 / 3 & r* = 50. 

In order to numerically perform the computation, we need to specify the exact expressions 
of the functions H(t) & /3(r) that appear in the expressions of section 2 9 

H(t) = /i(r) e 4Gl+4G2+2G3 & P(t) = J h{r) e 2G3 . (4.1) 

9 9 

As can be seen from both figures 2 & 3 the EE can either be positive or negative depending 
on the length of the strip. When the EE is positive the solution is a disconnected surface and 
when it is negative a connected. The transition between these two phases is characteristic 
in confining theories and first appeared in the literature in [14] (see also [15]). 

The connected solution exists only until a maximum value for the length of the strip 
(the point where L'(t) = 0 & L"{r) < 0). As can be seen from the plots of both figures 2 
& 3, there are two possible branches for the connected solution. The upper branch is an 
unstable solution while the lower is a stable one. The double valueness corresponds to a 
first order phase transition between the connected and the disconnected solution, at the 
critical point L cr a (it is the point at which the EE is zero). In [14] the authors claimed that 

8 See more about the effect of such a choice for the UV cutoff on the IR physical quantities in [21]. EE 
is even more less affected, since in (2.9) we are subtracting the contribution from the UV 

9 Note that the dilaton has disappeared from H(r) since we are working in the string frame. 
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Figure 2. On the left part of the figure it is the plot of the length of the strip as a function of 
To, which is the minimal value of r in the bulk of the geometry. In the right part of the figure it 
is the plot of the EE as a function of L. The black curve is for zero flavors and reproduces the 
computation of [14]. The red(blue) curve is for Nf = 0.1 (IV/ = 0.5). 

every change in the trend of the function L(tq ) is a sign for a confinement/deconfinement 
phase transition. In other words, every peak of the function L(to) corresponds to a phase 
transition. We will support this claim with our computations both in the massless & 
massive KS background. 

In figure 2 we plot the length of the strip and the EE as a function of tq as well as the 
EE as a function of L. As can be seen from the left plot of this figure, the critical length 
L C rit decreases as the number of flavors increases. Naively, from a weak field theory point 
of view, that is what we should expect since increasing the number of flavors will in turn 
increase screening. From another point of view increasing the number of flavors, increases 
the beta function and decreases the QCD scale A qcd- Since we expect the location of 
the transition (between the connected to the disconnected configuration) to happen at 
Lcrit ~ critical length should increase with the number of flavors. This was 

observed in [24] for the behavior of the critical length in the computation of the EE that 
was performed in the MN background with massless dynamical flavors [25]. The reason 
behind introducing an alternative way of fixing the constant h o (Case II) was to show that 
it is possible even in the massless KS background to reproduce such a behavior for the 
critical length. 

Before analyzing the behavior of the EE when we fix ho along the lines of Case II, we 
would like to point out a couple of things about the calculation of the EE in the (massless) 
MN vs. KS backgrounds and the importance of the value of ho for the appearance of a 
phase transition. Comparing the calculations for the EE between the KS and the MN 
backgrounds (without flavors), performed in [14] and [24] respectively, it is clear that 
while in the KS case there is a phase transition in the MN case there is not, since L is a 
monotonically increasing function. In [24] in order to reproduce confining behavior they 
were forced to introduce, by hand, a UV cut-off scale. On the other hand performing the 
computation of the EE in the KS background it is clear that the choice of fixing ho is very 
important. Had they chosen to fix it in a way that the warp factor is constant at infinity, 
the phase transition in the EE immediately disappears. 

Combining these two observations we conclude that if we want to fix ho along Case 
II and still be able to reproduce confinement, we should exploit the presence of the flavor 
driven UV cut-off scale. In this way confinement appears close to the UV cut-off scale, 
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in a way that seems artificial. There are ways to circumvent this by performing a UV 
completion of the background, along the lines that are extensively discussed in [15]. 



40 


30 


20 


10 


Figure 3. This is the analogue of figure 2, when we fix ho as it is described in Case II. The red(blue) 
curve is for Nf = 0.1 (Nf = 0.5) 

In figure 3 we plot the length of the strip and the EE as a function of to as well as 
the EE as a function of L, but now ho is fixed as it is described in Case II. In this way, 
contrary to what it is described in figure 2, the critical length increases and the glueball 
mass decreases with the number of flavors. The glueball scale is set by hand through the 
choice of r*. 

In the following paragraph we will study the EE in the massive KS background and 
focus the analysis on quantum phase transitions. For this reason we will not consider Case 
II and all the calculations will be performed in Case I. 

5 Entanglement entropy Sc massive flavors 

Moving one step forward from the setup we examined in the previous section, we will 
consider the strip in the Klebanov-Strassler with massive dynamical flavors and calculate 
its EE. The extra scale in the system, represented by the position that the source for the 
flavor branes gets activated, will introduce a second peak in the function L(t o) that will 
correspond to a another first order phase transition for the EE. The interplay/competition 
between the scales of the glueballs and the quarks will lead us to a critical point where 
one of the two phase transitions will become of the second order. At the point where the 
first order phase transition degenerates to a second it is possible to calculate the critical 
exponents and determine the universality class that the phase transition belongs. Even 
if we perform the calculation in the Heaviside approximation (3.13), the critical point we 
observe is genuine and not an artifact of the approximation. 

When dynamical flavors are present in the geometry, there is the notion of a screening 
length for the strip (in full analogy with the screening length in the Wilson loop compu¬ 
tation, e.g. [26],[20]). This is the length at which the connected surface would break by 
producing a pair of two cigars that end on a some of the smeared flavor branes. The value 
for that length would be determined by the condition that the EE for the connected surface 
equals the EE of two cigars ending on a flavor brane. Due to the smearing of the flavors 
the decay of the EE into such a channel is 1/IV C suppressed, and not Nf/N c as it would be 
for the case of the parallel localized flavor. 
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Since in the previous section we thoroughly examined the dependence of the EE to 
the number of flavors in this section we will focus on the effect of the quark mass. For this 
reason we fix the number of flavors to the value Nf = 1 and plot the function L(tq, T q = /x) 
(and the EE) for different values of the quark mass. When the two scales, namely the 
glueball and the quark scale, are noticeably far from each other (green vs. red plot in 
figure 4) the triangle (visualizing the first order phase transition) that corresponds to the 
quark scale is in the physical branch of the EE. Indeed in figure 4 the red triangle is in the 
unphysical region (for EE < 0), while the green triangle in the physical. 

We should mention that by comparing the plots of the EE with respect to the length of 
the stripe, for different number of flavors, we can see that the value of jx that corresponds 
to the critical value for the length increases when you increase the number of flavors. 


L.1CT 4 ' 3 



Figure 4. On the left part of the figure it is the plot of the length of the strip as a function of 
To for two different radii of the spherical cavity. The red (green) curve is for p, = 3.45 (/x = 3.65). 
In the right part of the figure it is the plot of the EE as a function of L. In order for the physical 
branch of the EE to contain two phase transitions (the lower branch of the second phase transition 
that corresponds to the disconnected solution is for EE = 0) the second peak should me noticeably 
lower than the first one, green vs. red (see also the analysis with a variety of examples in the 
appendix B.4 of [15]). Notice that as we decrease /x the spike related to the quark mass approaches 
and finally overtakes the one related to the glueball scale. The calculation is for Nf = 1 . 

Keep decreasing the distance between the two energy scales we come to a point where 
one of the two peaks (the one that corresponds to the glueballs) is dissolved and we are 
left with just one peak (see the left part of figure 5). This peak corresponds to the energy 
scale set by the mass of the quarks. The critical value of /x where the first peak disappears 
is a transition point between a first and a second order phase transition. The calculation 
of the corresponding critical exponents will put us in the position to make a conjecture 
about the existence of a universality class for the transitions of the EE 10 . 

In order to study and characterize the phase transition we will focus our analysis on 
the behavior of the function L(to,/u). Plotting this function we are looking for the critical 
value of /x where not only the first but also the second derivative of the function L(ro,/x) 
with respect to tq vanishes. Near the critical point the function L(tq,/x) will be well 

1(, The calculation of the critical exponents is performed using the Heaviside approximation (3.13) for 
the flavor distribution function. Because of this, as can be seen from the left plot of figure 4, the peak 
that corresponds to the quark scale is more acute than the peak of the glueball scale. Since the peak that 
disappears, and for which we calculate the critical exponents, is the smooth one, the use of the heaviside 
approximation has no effect on the values of the exponents and the determination of the universality class. 
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Figure 5. On the left part of the figure it is the plot of the length of the strip for three different 
radii of the spherical cavity. In the blue curve we have used the critical value of //, while the red 
(green) curve is for values of /i above (below) this critical value. The values we have used are 
/r = 3.2 & [x = 2.5. In the right part of the figure it is the plot of the EE as a function of L. As 
can be seen from this plot the second triangle disappears exactly at the critical value of /r and if 
we keep on decreasing the critical value we are left with just one phase transition. The calculation 
presented in this plot is for Nf = 1 . 


approximated by the following expression 


L - L c 


Oo - r c f . 


(5.1) 


Performing the calculation after fixing the number of flavors to Nf = 1 we get the following 
values for the critical length, fx & r 

L c » 10.62 t c » 2.59 & n c « 2.92. (5.2) 


From the behavior of the critical curve it is possible to extract the value for the critical 
exponents (for their definition see [27]). Since along the critical value /x = fx c the behavior 
of the critical curve is L — L c ~ (to — r c ) 5 , the value for this critical exponent is <5 = 3. 
In order to specify the critical exponent that is called (3, we need to fit the numerical data 
to the following curve 


to - t c 


t c 



(5.3) 


for fx < [x c ft t L = L c . As can be seen from figure 6, a non-linear fit (the solid blue curve) 
to the numerical data of the parametric plot of the minimal radial position in the bulk as 
a function of the quark mass, provides the critical value /3 = 1/3 with a high accuracy. In 
order to explore the possible dependence of this critical exponent to the number of flavors 
Nf, we performed the same computation for three different values of Nf, namely 1/2, 1 
& 3/2 11 . Our findings to high accuracy indicate that the value (3 = 1/3 is universal and 
independent of the number of flavors. Although the critical value for 5 is the classical one 
(in agreement with the classical mean field theory value) the critical value for (3 is not. So 
far in the literature there are examples of calculations of critical exponents (in a variety of 
gravity backgrounds) that either agree (see e.g [9]) or disagree (see e.g. [28, 29]) with the 
expectations that originate from the classical mean field theory value. Thus, overall the 


11 In figure 6 we present the numerical fit for Nf = 1 (Nf 


1 /2) corresponding to the left (right) plot. 
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behaviour of the strip width around the critical point takes the following form 



(5.4) 


Notice the difference between (5.4) and equation (5.1) of [9]. In the latter the authors 
studied the heavy quark antiquark potential (expectation value of two Wilson lines) in 
the massive KS background. They found the phase transition in the universality class of 
the Van der Waals system with mean field theory critical exponents. They argued that 
this is related to the theory of singularity of families of functions, or Catastrophe theory, 
and that this universality reflects the fact that they were analyzing a classical object: a 
macroscopic string. In our case the basic object is still classical, a minimal surface, but 
the system belongs in a different universality class, probably because we are dealing with a 
higher dimensional surface. A final comment is in order. Let us stress that the first order 
phase transition of [9, 26] which is related to the existence of dynamical quarks is present 
only for values of // smaller than a critical value. When /r is greater than this critical value 
their phase transition ceases to exist. In our case the quark related phase transition is 
present for all values of [i becoming dominant for ^ < /x cr . However, one should always 
keep in mind that we are considering different observables. We focus on the EE while the 
aforementioned studies focus on the heavy quark potential. 

Two of the critical exponents, through the use of the following scaling relations 


a + 2/3 + 7 = 2 , 7 = (3(5 — 1) , 2 — a = ud , 7 = u(2 — r]) , (5-5) 

are sufficient to completely determine the universality class of the system, in case we know 
its dimension d. This dimension is related to the computation of the correlation function 
and since this calculation is beyond the scope of the paper, we will restrict on the critical 
exponents a $z 7 . Using (5.5) and the already known values for 6 (3 we have 

a = 7 = ^, P = \ & 5 = 3. (5.6) 




Figure 6. In the two plots of that figure we represent (with the solid blue curve) a non-linear fit 
to the data of the numerical computations for the minimal radial position in the bulk as a function 
of the quark mass. On the left (right) plot the data are for Nf = 1 (Nf = 1/2). In both cases the 
value of the critical exponent is universal, independent of the value of Nf and different from the 
classical mean field theory value, namely (3 = 1/3. 
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Finally, let us comment on the behaviour of the EE as a function of the width of 
the strip L around the critical point. As discussed above at fi = fi c the first order phase 
transition becomes a second order one with tq playing the role of the order parameter. If 
we wanted to make a comparison with a thermodynamic system, the EE corresponds to the 
Gibbs free energy, the length of the strip corresponds to the pressure, the order parameter 
ro corresponds to the density while the quark mass /i corresponds to the temperature. As 
in the the quark antiquark potential one can use (2.5) and (2.9) to prove the following 
relations 

k = jm . (5.7) 

dL 2 dL 2 ±y/Hfa)L'(Tb) 

Near the critical point the function L(tq) is well approximated by (5.1). Substituting the 
derivative of (5.1) in the second relation of (5.7) and replacing to — r c in terms of L — L c 
we can integrate to get the following expression 


S - S c 


^/HW) {l 




(5.8) 


where D is a constant that depends on L, H and their derivatives at the critical point r c . 


6 Conclusions &; outlook 

In this paper we studied the EE of the confining KS background, when it is coupled to a 
large number of dynamical flavors (massive & massless) in the Veneziano limit. Dividing 
one of the spatial directions into a line segment of length L and its complement the EE 
between the two subspaces admits a first order phase transition, when we vary the width 
of the strip L. 

In section 4 we analyze the behavior of both the length of the strip and the EE as we 
vary the number of massless flavors. Depending on the way we fix the integration constant 
of the warp factor, that is related to the glueball scale, we find two different quantitative 
behaviors. Requiring the warp factor to vanish at some finite UV cut-off (in analogy with 
[7]), the glueball mass increases with the number of flavors. The opposite behavior for the 
glueball mass occurs, when we impose on the warp factor to be the same at the origin for 
any value of Nf. 

In section 5 the mass of the flavor branes introduce another scale in the background. 
In this case the plot of the length of the strip presents a second peak at the position 
that the flavor branes gets activated and the EE presents a second phase transition. The 
competition between the glueball and the quark scales will lead to a critical point where 
the phase transition which was originally associated with the glueball scale will become a 
second order one. The critical exponents we calculate are independent of the number of 
flavors and different from the mean filed theory expectations. This is in contradistinction 
with similar calculations of the quark-antiquark potential in the same background in [9]. 
On general grounds, since we consider the effect of the fundamental loops in the Feynman 
diagrams of the gauge theory, such a deviation from the mean field theory expectations 
for the critical exponents should not be unexpected. Even in the quenched approximation 
quantum effects are sufficiently strong for that, see e.g. [29]. 

Finally, let us briefly comment on the physical picture behind our findings. When the mass 
of the quarks is much larger that the lightest glueball scale, i.e. when /r >> fi cr , the first 
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order transition related to the existence of the quarks is happening closer to the UV cut-off 
and affects little the IR physics, that is the QCD scale A qcd or the lightest glueball mass. 
But as we decrease the quark mass the effect of the quarks becomes significant since we 
can no longer integrate them out. Indeed, as can be seen from Figures 4 and 5 when we 
decrease /i the critical length L cr increases and as a result the QCD scale A qcd 1/L 

cr 

gets smaller. This qualitative behaviour can be understood if we consider the two extreme 
limits. When /r fjL cr the quarks are very heavy and they decouple from the dynamics 
of the theory. Effectively this corresponds to setting Nf = 0. On the other hand, when 
fj, ~ n cr the presence of a large number of quarks increases the /3-function which results in 
a smaller QCD scale. 

We will close this section by describing a couple of ideas for future studies. In [15] they 
presented a detailed study of the EE for several confining theories, with and without flavor. 
Some of the gravity models they examine, present multiple phase transitions in the EE, 
due to the presence of various scales (e.g. quarks with different masses). One could try to 
study the behavior of the gravity system when the scales (masses) become comparable, by 
calculating the critical exponents. It would be interesting to check whether those systems 
belong to the same universality class with the one we describe in this paper. Since both 
peaks of the length of the strip will correspond to quark masses, another calculation would 
be to examine the possible dependence of the critical exponents on the number of flavors 
of the different spherical cavities. 

In [22] the construction of a gravity background with backreacted flavor branes of dif¬ 
ferent masses in the AdS 5 x S ' 5 was presented. The flavor branes produce multiple spherical 
cavities with radii that correspond to the quark masses. Generalizing that construction to 
the KS background will imply a confining model with multiple scales and multiple phase 
transitions. It would be interesting to study the interplay of the different scales in the EE 
and the values of the critical exponents. 
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